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^ : Abstract 
O 

The proton-neutron relativistic quasiparticle random phase approximation (PN-RQRPA) is for- 
mulated in the canonical single-nucleon basis of the relativistic Hartree-Bogoliubov (RHB) model, 

O , 

^ , for an effective Lagrangian characterized by density-dependent meson-nucleon couplings. The 



model includes both the T = 1 and T = pairing channels. Pair configurations formed from the 



^ ' fully or partially occupied states of positive energy in the Fermi sea, and the empty negative-energy 

H : 

states from the Dirac sea, are included in PN-RQRPA configuration space. The model is applied to 
the analysis of charge-exchange modes: isobaric analog resonances and Gamow- Teller resonances. 

PACS numbers: 21.30.Fe, 21.60.Jz, 24.30.Cz, 25.40.Kv 
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I. INTRODUCTION 



A consistent and unified treatment of mean-field and pairing correlations is crucial for a 
quantitative analysis of ground-state properties and multipole response of unstable, weakly 
bound nuclei far from the line of /5-stability. In Ref. |l| we have formulated the relativistic 
quasiparticle random phase approximation (RQRPA) in the canonical single-nucleon basis of 
the relativistic Hartree-Bogoliubov (RHB) model. The RHB model presents the relativistic 
extension of the Hartree-Fock-Bogoliubov framework, and it provides a unified description 
of particle-hole [ph) and particle-particle {pp) correlations. In this framework the ground 
state of a nucleus can be written either in the quasiparticle basis as a product of indepen- 
dent quasi-particle states, or in the canonical basis as a highly correlated BCS-state. By 
definition, the canonical basis diagonalizes the density matrix and it is always localized. It 
describes both the bound states and the positive-energy single-particle continuum. The for- 
mulation of the RQRPA in the canonical basis is particularly convenient because, in order to 
describe transitions to low-lying excited states in weakly bound nuclei, the two-quasiparticle 
configuration space must include states with both nucleons in the discrete bound levels, 
states with one nucleon in a bound level and one nucleon in the continuum, and also states 
with both nucleons in the continuum. 

The .e,ativ.tic QRPA of Ref. [l] . fully .lf-co,.i.ten. Fo. the i,.e.aetio„ i„ the pattHe- 
hole channel effective Lagrangians with nonlinear meson self-interactions are used, and pair- 
ing correlations are described by the pairing part of the finite range Gogny interaction. 
Both in the ph and pp channels, the same interactions are used in the RHB equations that 
determine the canonical quasiparticle basis, and in the matrix equations of the RQRPA. 
The RQRPA configuration space includes also the Dirac sea of negative energy states. The 
RHB+RQRPA approach has been tested in the example of multipole excitations of neutron 
rich oxygen isotopes, and the model has been employed in the analysis of the evolution of 
the low-lying isovector dipole strength in Sn isotopes and N=82 isotones. 

Relativistic mean-field and RPA calculations based on effective Lagrangians with non- 
linear meson self-interactions present not only a number of technical problems, but also 
the description of finite nuclei obtained with these effective interactions is not satisfactory, 
especially for isovector properties. As an alternative, in a number of recent applications 
models with density-dependent meson-nucleon vertex functions have been used. Although 
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the two classes of models are essentially based on the same microscopic structure, i.e. on 
density-dependent interactions, the latter can be more directly related to the underlying 
microscopic nuclear interactions. Several recent analyses have shown that relativistic effec- 
tive interactions with explicit density dependence of the meson-nucleon couplings provide 



an improved descri 
stability. In Ref. j2 



)tion of asymmetric nuclear matter, neutron matter and nuclei far from 
we have extended the relativistic Hartree-Bogoliubov (RHB) model 
to include density dependent meson-nucleon couplings. The effective Lagrangian is charac- 
terized by a phenomenological density dependence of the cr, uo and p meson-nucleon vertex 
functions, adjusted to properties of nuclear matter and finite nuclei. It has been shown that, 
in comparison with standard RMF effective interactions with nonlinear meson-exchange 
terms, the new density-dependent meson-nucleon force DD-MEl significantly improves the 
description of asymmetric nuclear matter and of ground-state properties oi N ^ Z nuclei. 
This is, of course, very important for the extension of RMF-based models to exotic nuclei 
far from /3-stability, and for applications in the field of nuclear astrophysics. 

In Ref. j^l the matrix equations of the relativistic random-phase approximation (RRPA) 
have been derived for an effective Lagrangian characterized by density-dependent meson- 
nucleon couplings. The explicit density dependence introduces rearrangement terms in the 
residual two-body interaction, and their contribution is essential for a quantitative descrip- 
tion of excited states. Illustrative calculations of the isoscalar monopole, isovector dipole 
and isoscalar quadrupole response Of208pb, were performed in the fully self-consistent RRPA 
framework based on effective interactions with a phenomenological density dependence ad- 
justed to nuclear matter and ground-state properties of spherical nuclei. 

In this work we extend the model developed in Ref. ^|, and formulate a proton- neutron 
relativistic QRPA in the canonical single-nucleon basis of the RHB model. The proton- 
neutron RQRPA, with the density-dependent effective interaction DD-MEl, will be em- 
ployed in the analysis of nuclear spin and isospin excitations. 

Collective spin and isospin excitations in atomic nuclei have been the subject of many 
experimental and theoretical studies (for an extensive review see Ref. ^). Nucleons with 
spin up and spin down can oscillate either in phase (spin scalar S=0 mode) or out of phase 
(spin vector S=l mode). The spin vector, or spin-flip excitations can be of isoscalar (S=l, 
T=0) or isovector (S=l, T=l) nature. These collective modes provide direct information 
on the spin and spin-isospin dependence of the effective nuclear interaction. 
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Especially interesting is the collective spin-isospin oscillation with the excess neutrons 
coherently changing the direction of their spins and isospins without changing their orbital 
motion - the Gamow- Teller resonance (GTR) J'^ = 1+. This collective mode was predicted 

n 

already in 1963 |5|, but it was only in 1975 that the first experimental indications of the 
GT resonance were observed in {p, n) charge-exchange reactions at intermediate energies . 
The detailed knowledge of GT strength distributions is essential for the understanding of 
nuclear reactions in the process of nucleosynthesis. In particular, the low-lying GT strength 
is directly related to /3-decay rates, as well as to the electron-capture process leading to the 
stellar collapse and supernovae explosion. In many nuclei all over the periodic table the 
GT strength distribution, when measured in the excitation energy region where the most 
pronounced GT peaks occur, is quenched by more than 20% when compared to a model 
independent sum rule. Two physically different mechanisms had been suggested as a possible 
explanation of the quenching of the total GTR strength: (i) nuclear configuration mixing 
- the high-lying 2p — 2h states mix with the Ip — Ih GT states and shift the GT strength 
to high-energy region far beyond the resonance jz, 0]; (ii) the high-energy A-isobar - 
nucleon-hole configurations (A — h) couple to the GT mode and remove part of the strength 



from the low-lying excitation spectrum 



ll| . Recent {p, n) scattering experiments have 



shown, however, that only a small fraction of the GT quenching originates from A — h 
transitions jl^ . 

The GT spin-flip isovector excitations and the related /3-decays have been theoretically 
investigated by employing: (i) the shell-model approach, (ii) the non-relativistic proton- 
neutron quasiparticle random phase approximation, and (iii) the relativistic random phase 
approximation in closed-shell nuclei. Experimental data on charge-exchange excitations in 
light and medium-mass nuclei are very successfully reproduced by large-scale shell-model 
calculations. However, as the number of valence nucleons increases, the dimension of shell- 
model configuration space becomes far too large for practical calculations. Shell-model 
calculations have been recently extended to the region of p/-shell nuclei with A = 45 — 
65. Experimental Gamow- Teller strength distributions and nuclear /3-decay half-lives have 
been reproduced The GT response in medium- mass nuclei has also been successfully 
described by shell-model Monte Carlo calculations (SMMC) jl4| . 

The proton - neutron quasiparticle random phase approximation (PN-QRPA) can be 
used to describe nuclei in mass regions that are presently beyond the reach of the most 
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sophisticated shell-model calculations. The pioneering work of Halbleib and Sorensen Q| 
was based on simple separable forces. Further developments employed zero-range interac- 
tions with BCS-type pairing. Most studies of GT excitations have been based on Skyrme 

It has also been shown that the inclusion of particle-particle 
correlations in the QRPA residual interaction has an important effect in the calculation of 
/?-decay rates H H B, Q, E , 26 1 . 

Experimental data on unstable nuclei close to N=Z line, that became available in recent 
studies with radioactive beams, have also renewed the interest in the role of proton-neutron 
pairing. In medium-heavy and heavy nuclei closer to the line of /3-stability, the proton and 
neutron Fermi levels are located in different major shells, and therefore the contribution of 
proton-neutron pairing correlations to the ground-state binding is usually neglected. The 
low-lying GT strength responsible for /3-decay, however, involves proton-neutron particle- 
particle correlations and is sensitive to the T = pairing interaction. In Ref. /3-decay 
rates for spherical neutron-rich r-process waiting-point nuclei have been calculated within 
a fully self-consistent QRPA, formulated in the Hartree-Fock-Bogoliubov canonical single- 
particle basis. It has been shown that the proton-neutron particle-particle interaction has a 
pronounced effect on the calculated half-lives. 

Surprisingly little work has been reported on the description of charge-exchange excita- 
tions in finite nuclei in the framework of relativistic mean-field models. The first relativistic 
RPA calculations of isobaric analog resonances (lAR) and Gamow- Teller resonances have 

n n 

been performed only recently in Refs. |22l, U^- This analysis was, however, restricted to 



doubly closed-shell nuclei. A rather small configuration space was used and, furthermore, 
co„flgu.at,o„= that Mude empty states tor. the negative-energy Dnac sea were neglected. 
In a very recent work f29] it has been shown, on the other hand, that the inclusion of 
negative-energy Dirac sea states has a pronounced effect on the calculated Gamow- Teller 
sum rule. A more complete relativistic RPA calculation of GT resonances in doubly closed- 
shell nuclei has been reported in Ref. j^^. The ground states of ^'^Ca, ^°Zr and ^osp^ -were 
calculated in the relativistic mean-field (RMF) model with the NL3 effective interaction. 
The spin-isospin correlations in the RRPA calculation of the GT response functions were in- 
duced by the isovector mesons vr and p. In addition to the standard non-linear NL3 effective 
interaction j^| with the vector rho-nucleon coupling, the effective Lagrangian included the 
pseudo-vector pion-nucleon interaction. Although in the relativistic mean-field description 
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of the nuclear ground state the direct one-pion contribution vanishes at the Hartree level 
because of parity conservation, the pion nevertheless plays an important role for excitations 
that involve spin degrees of freedom. Since it has a relatively small mass, the pion mediates 
the effective nuclear interaction over large distances. The p — h residual interaction with 
p- and vr-meson exchange has been extensively used in non-relativistic RPA calculations of 



charge-exchange excitations 



32 



33|. Because of the derivative type of the pion-nucleon cou- 



pling, it is also necessary to include a zero-range Landau-Migdal term that accounts for the 
contact part of the nucleon-nucleon interaction. The analysis of Ref. 30| has shown that the 
RRPA calculation with the NL3 effective interaction, the pseudo-vector pion-nucleon cou- 
pling (m^ = 138 MeV and /^/47r = 0.08), and the Landau-Migdal force with the strength 
parameter Qq = 0.6, reproduces the experimental excitation energies of the main components 
of the GT resonances in ^^Ca, ^°Zr and ^°®Pb. 

In Ref. we have calculated the GTR and lAR for a sequence of even-even Sn target 
nuclei by using the framework of the relativistic Hartree-Bogoliubov model plus proton- 
neutron quasiparticle random-phase approximation. The calculation reproduces the exper- 
imental data on ground-state properties, as well as the excitation energies of the isovector 
excitations. It has been shown that the isotopic dependence of the energy spacings between 
the GTR and lAR provides direct information on the evolution of neutron skin-thickness 
along the Sn isotopic chain. A new method has been suggested for determining the difference 
between the radii of the neutron and proton density distributions along an isotopic chain, 
based on measurement of the excitation energies of the GTR relative to the lAR. 

This work presents a much more detailed analysis of charge-exchange modes, both for 
doubly closed-shell and for open-shell nuclei. In Sec. [H] we introduce the formalism. The 
matrix equations of the proton-neutron relativistic QRPA are formulated in the canonical 
basis of the relativistic Hartree-Bogoliubov (RHB) model for spherical nuclei. In Sec. Illll the 
RHB+RQRPA model is employed in illustrative calculations of charge-exchange collective 
modes: isobaric analog resonances and Gamow- Teller resonances. Section IIVI contains the 
conclusions. 
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II. PROTON-NEUTRON RELATIVISTIC QUASIPARTICLE RANDOM PHASE 
APPROXIMATION 

In Ref. Q] we have derived the relativistic quasiparticle random phase approximation 
(RQRPA) from the time- dependent relativistic Hartree-Bogohubov (RHB) model in the 
limit of small amplitude oscillations. The full RQRPA equations are rather complicated 
and it is considerably simpler to solve these equations in the canonical basis, in which the 
relativistic Hartree-Bogoliubov wave functions can be expressed in the form of BCS-like 
wave functions. In this case one needs only the matrix elements V^Vk'a of residual ph- 
interaction, and the matrix elements V^J!,-^y of the pairing pp-interaction, as well as certain 
combinations of the occupation factors f The details of the RHB model, of the RQRPA 
formalism, and the technicalities of the solution of the RQRPA equations in the canonical 
basis, are given in Ref. jlj]. In this section we only collect those expressions which are 
specific for the proton-neutron relativistic quasiparticle random phase approximation (PN- 
RQRPA). The model can be considered as a relativistic extension of the fully-consistent 
proton-neutron QRPA formulated in Ref. and applied in an analysis of /3-decay rates 
of r-process nuclei. 

We consider transitions between the 0"^ ground state of a spherical even-even parent 
nucleus and the J'^ excited state of the corresponding odd-odd daughter nucleus. These 
transitions are induced by a charge-exchange operator T'^^^ . Taking into account the ro- 
tational invariance of the nuclear system, the quasiparticle pairs can be coupled to good 
angular momentum and the matrix equations of the PN-RQRPA read 



J ^"j I 1 yAj I I -1 j 



The matrices A and B are defined in the canonical basis (35i] 



(1) 



+ {UpVnUp'Vn' + fpM„fp/M„/) V^'^,'^^, + {UpUnUp'Un' + VpVnVp'Vn') V^^^,^, 
Bin,p'n' = {-^y^'^^-'^'^UpVnVp'Un' +VpUnUp,Vn')V^^^^, 

- {UpUnVp'Vn' + VpVnUp'Un') V^^^,^, . (2) 

Here p, p', and n, n' denote proton and neutron quasiparticle canonical states, respectively. 



V^^ is the proton-neutron particle-hole residual interaction, and V^^ is the corresponding 
particle-particle interaction. The canonical basis diagonalizes the density matrix and the 
occupation amplitudes Vp^n are the corresponding eigenvalues. The canonical basis, however, 
does not diagonalize the Dirac single-nucleon mean- field Hamiltonian hf) and the pairing 
field A, and therefore the off-diagonal matrix elements H^\, and 11^^, appear in Eq. (j2)): 

Hll, = {u^u^> - v^v^^)h^^: - {u^v^^ + v^u^:)A^^> , (3) 

For each energy Ex, X^^ and Y^^ in Eq. Q denote the corresponding forward- and 
backward-going QRPA amplitudes, respectively. The total strength for the transition be- 
tween the ground state of the even-even (N,Z) nucleus and the excited state of the odd-odd 
(N+1,Z-1) or (N-1,Z+1) nucleus, induced by the operator T"^^^ , reads 

2 



pn 



(4) 



The discrete strength distribution is folded by the Lorentzian function 



In the illustrative calculations included in the following sections, the choice for the width of 
the Lorentzian function is 1 MeV. 

The spin-isospin dependent interaction terms are generated by the p- and vr-meson ex- 
change. Because of parity conservation, the one-pion direct contribution vanishes in the 
mean-field calculation of a nuclear ground state. Its inclusion is important, however, in cal- 
culations of excitations that involve spin and isospin degrees of freedom. The particle-hole 
residual interaction in the PN-RQRPA is derived from the Lagrangian density 

/:r+p = -gpi'^P^T^ - —^Pi^rd.nf^ . (6) 

Vectors in isospin space are denoted by arrows, and boldface symbols will indicate vectors 
in ordinary three-dimensional space. 

The coupling between the p-meson and the nucleon is assumed to be a vertex function of 



the vector density = ^jfij^, with j^j = ip'j^ip. In Ref. Pi it has been shown that the explicit 
density dependence of the meson-nucleon couplings introduces additional rearrangement 
terms in the residual two-body interaction of the RRPA, and that their contribution is 
essential for a quantitative description of excited states. However, since the rearrangement 
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terms include the corresponding isoscalar ground-state densities, it is easy to see that they 
are absent in the charge exchange channel, and the residual two-body interaction reads 



V{ri,r2) = TiT2{[3Y)i{htJi)2gp{pv{ri))gp{p,{r2))Dp{ri,r2) 

2 

I rir2(SiVi)(S2V2)/^.(n,r2) . (7) 



Dp(^T,) denotes the meson propagator 

X g-»n-p(,r)kl-T-2 



Dpij,) = -. . 7— , (8) 

An \ri — r2\ 



and 

^ cr 



(9) 



a- 

For the p- meson coupling we adopt the functional form used in the DD-MEl density- 
dependent effective interaction 

9pipv) = gpipsat)exp[-ap{x - 1)] , (10) 

where x = Pv/psat, and psat denotes the saturation vector nucleon density in symmetric 
nuclear matter. For the pseudovector pion-nucleon coupling we use the standard values 

= 138.0 MeV — = 0.08 . (11) 

The derivative type of the pion-nucleon coupling necessitates the inclusion of the zero-range 
Landau-Migdal term, which accounts for the contact part of the nucleon-nucleon interaction 

Vs. = g'(^] rir2Si • S25(ri - r2) , (12) 

with the parameter g' ^ 0.6 adjusted to reproduce experimental data on the GTR excitation 
energies. 

With respect to the RHB calculation of the ground state of an even-even nucleus, the 
charge-exchange channel includes the additional one-pion exchange contribution. The PN- 
RQRPA model is fully consistent: the same interactions, both in the particle- hole and 
particle-particle channels, are used in the RHB equation that determines the canonical 
quasiparticle basis, and in the PN-RQRPA equation (Q). In both channels the same strength 
parameters of the interactions are used in the RHB and RQRPA calculations. 
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The two-quasiparticle configuration space includes states with both nucleons in the dis- 
crete bound levels, states with one nucleon in the bound levels and one nucleon in the 
continuum, and also states with both nucleons in the continuum. In addition to the con- 
figurations built from two-quasiparticle states of positive energy, the RQRPA configuration 
space contains pair-configurations formed from the fully or partially occupied states of pos- 
itive energy and the empty negative-energy states from the Dirac sea. As will be shown in 
the next section, the inclusion of configurations built from occupied positive-energy states 
and empty negative-energy states is essential for the consistency of the model. 

Nuclear properties calculated with the RHB plus RQRPA model depend on the choice 
of the effective RMF Lagrangian in the p/i-channel, as well as on the treatment of pairing 
correlations. In this work, for the RHB calculation of ground states and in the p/i-channel 
of the residual interaction, we use the density-dependent effective interaction DD-MEl. 
In Ref. Q the RHB model with the density- dependent interaction DD-MEl in the ph- 
channel, and with the finite range Gogny interaction DIS in the pp-channel, has been tested 
in the analysis of the equations of state for symmetric and asymmetric nuclear matter, 
and of ground-state properties of the Sn and Pb isotopic chains. It has been shown that, 
as compared to standard non-linear relativistic mean-field effective forces, the interaction 
DD-MEl has better isovector properties and therefore provides an improved description of 
asymmetric nuclear matter, neutron matter and nuclei far from stability. 

In the pp-channel of the RHB model we have used a phenomenological pairing interaction, 
the pairing part of the Gogny force. 



ypp(l,2) = e"""!-"^)/'^^)' {W, + BiP" - HiP^ - MiP^P^), (13) 

«=1,2 

with the set DIS 



36l | for the parameters /ij, Wj, -Bj, Hi and Mj {i = 1,2). This force has been 
very carefully adjusted to the pairing properties of finite nuclei all over the periodic table. 
In particular, the basic advantage of the Gogny force is the finite range, which automatically 
guarantees a proper cut-off in momentum space. In the present analysis we will also use 
the Gogny interaction in the T = 1 pp-channel of the PN-RQRPA. For the T = proton- 
neutron pairing interaction in open shell nuclei we employ a similar interaction: a short-range 
repulsive Gaussian combined with a weaker longer-range attractive Gaussian: 

^i2 = -K)E^7,e-'--/^'n5=i,T=o , (14) 
i=i 

10 



where n5=i,T=o projects onto states with S = 1 and T = 0. This interaction was used 
in the non-relativistic QRPA calculation j3l of /9-decay rates for spherical neutron-rich r- 



process waiting-point nuclei. As it was done in Ref. j26|], we take the ranges /xi=1.2fm and 
/i2=0.7fm of the two Gaussians from the Gogny interaction ()13|) . and choose the relative 
strengths gi = 1 and g2 = —2 so that the force is repulsive at small distances. The only 
remaining free parameter is Vq, the overall strength. 

III. CHARGE-EXCHANGE COLLECTIVE MODES OF EXCITATIONS 
A. The isobaric analog resonance 

The isobaric analog resonance (lAR) presents the simplest charge-exchange excitation 
mode, detected already forty years ago in experiments on low-energy proton elastic scat- 
tering from heavy nuclei. The observed resonances were at energies consistent with the 
interpretation that they were isobaric analog states in the compound nucleus. One of the 
important characteristics of an lAR is its narrow width. This is because it has the same 
isospin as the parent state, while the neighboring states have the isospin of the ground state 
of the daughter nucleus, i.e. they differ in isospin by one unit. This means that they will 
couple only weakly with the lAR. 

As a first test of our PN-RQRPA model, we calculate the lAR strength functions, i.e. 
the distribution of J'^ = 0~^. The Fermi transition operator reads 

Til = E^±- (15) 

i=l 

In Fig. ^we display the PN-RRPA response to the operator (jl5|) for ^^Ca, ^°Zr and ^°^Pb. 
The strength distributions are dominated by a single lAR peak, which corresponds to a 
coherent superposition of np — uh (or proton-neutron 2qp) excitations. The calculated 
excitation energies (evaluated with respect to the ground state of the parent nucleus) are 
compared with the _nding expartoental value, (thick a.ows) from (p,n) scattering 
data for ^^Ca j37|, ^°Zr |]J,|38|, and ^o^Pb [39]. The agreement between the PN-RRPA and 
experimental data is indeed very good. We have also verified that the calculated strength 
distributions exhaust the Fermi non-energy weighted sum rule 

Sf,-{F)=N-Z. (16) 

11 



The proton-neutron relativistic QRPA provides a natural framework for the description 
of spin and isospin excitations in open-shell nuclei. In Fig.|21we plot the calculated lAR exci- 
tation energies for the sequence of even-even Sn target nuclei with A = 108 — 132. The result 
of fully self-consistent RHB plus proton-neutron RQRPA calculations are shown in compar- 
ison with experimental data obtained in a systematic study of the (^He,t) charge-exchange 
reaction over the entire range of stable Sn isotopes The calculated values reproduce the 
empirical mass dependence of the lAR excitation energies, and we also notice a very good 
agreement with available experimental data. For the nuclei with A = 112 — 124 the largest 
difference between the theoretical and experimental lAR excitation energies is ~ 200 keV. 
For the Sn nuclei with A = 108 — 132, in Figs. 121 and H] we display the corresponding theoret- 
ical strength distribution functions. To illustrate the importance of a consistent treatment 
of pairing correlations, in addition to the response functions calculated with the fully self- 
consistent RHB plus proton-neutron RQRPA, we also display the strength functions obtained 
without including the proton-neutron residual pairing interaction. Proton-neutron pairing, 
of course, does not contribute in the RHB calculation of the ground states. The dashed 
curves are calculated by including only the p/i-channel of the RQRPA residual interaction. 
In both cases (solid and dashed curves) T = 1 neutron-neutron pairing has been included in 
the RHB calculations of the ground-states of Sn target nuclei. The results shown in Figs. El 
and m clearly illustrate the pronounced effect that T = 1 proton-neutron pairing correlations 
have on the calculated lAR. Without proton-neutron pairing the calculated lAR excitation 
energies for nuclei with A> 116 are simply too high. The attractive proton- neutron pairing 
interaction redistributes the strength between various QRPA components, and lowers the 
calculated excitation energy to the value which corresponds to the difference between the 
Coulomb energies of the target nucleus and the {Z + 1, N — 1) daughter nucleus. As one 
would expect, the effect gradually weakens as the closed shell at = 82 is approached. The 
effect of the proton-neutron residual pairing on the lAR strength distributions in lighter 
Sn isotopes with A < 114 is more dramatic. Without the contribution of the pp-channel 
to the proton-neutron residual interaction, the response function is fragmented and it does 
not display the experimentally observed single narrow resonance. Only with the inclusion 
of proton-neutron pairing a single collective state is obtained. 

The mechanism through which the PN- RQRPA builds the lAR in Sn isotopes is illustrated 
in Fig. El for ^"'^Sn, ^^'^Sn and ^^°Sn. The figure displays the discrete QRPA spectra for the 
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following three cases. In the panels on the left we plot the spectra calculated without 
any pairing interaction, either in calculation of the ground-state, or in the residual proton- 
neutron interaction. It is immediately clear why the three particular examples were selected. 
Without pairing, the (77/2 neutron orbital is fully occupied in the ground-state of ^°^Sn, the 
(^7/2 and the c/5/2 orbitals are fully occupied in the ground-state of ^^''Sn, and for ^^°Sn all 
the neutron orbitals in the shell = 50 — 82 are occupied, except /iii/2. For all three nuclei 
the response function to the Fermi operator p5|) displays a single peak at the excitation 
energy very close to the experimental position of the lAR. The distribution of neutron-to- 
proton RPA amplitudes is shown next to each major peak. For example, in the case of 
^°*Sn 99% of the strength results from the transition v (77/2 vr (77/2, whereas for ^^°Sn 
98% of the strength is distributed between the four orbitals occupied by the neutrons in 
the = 50 — 82 shell. In the next step (middle panels in Fig. we include the Gogny 
neutron-neutron pairing in the RHB calculation of the ground states, and as a result all 
the A^ = 50 — 82 neutron orbitals acquire finite occupation probabilities. The pp-channel, 
however, is not allowed to contribute to the QRPA matrix elements of the proton-neutron 
residual interaction. In this case the QRPA spectra become fragmented in ^°^Sn and ^^^Sn. 
Considerable strength is transferred to states which correspond to transitions from neutron 
orbitals with low occupation probabilities: S1/2 and /iii/2. For the most pronounced peaks we 
have included the distribution of neutron-to-proton QRPA amplitudes. The fragmentation 
does not occur in ^^°Sn where, as a result of ground-state pairing correlations, all the A^ = 
50 — 82 neutron orbitals have considerable occupation probabilities. In ^^"Sn we notice only 
a redistribution of QRPA amplitudes in favor of higher-lying orbitals, and consequently the 
repulsive particle-hole residual interaction shifts the lAR to higher energy. For the two 
nuclei from the lower half of the A^ = 50 — 82 neutron shell, the 'ph residual interaction 
fragments the QRPA spectra and effectively precludes the formation of the lAR. The panels 
on the right-hand side of Fig. El display the result obtained by fully-consistent PN-RQRPA 
calculations, including the T = 1 proton-neutron pairing channel: the strength has been 
collected in a single peak - the lAR. We notice that with respect to the case (a) without any 
pairing correlations, the lAR strength has been redistributed among the neutron-to-proton 
QRPA amplitudes, reflecting the occupation of higher-lying neutron orbitals in the parent 
nucleus. The calculated excitation energies, as we have shown in Fig. |21 are in very good 
agreement with the experimental data on the lAR in Sn nuclei Q|. The total effect is, of 
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course, simply a consequence of the fact that the pairing interaction is isospin invariant, and 
therefore commutes with the Fermi operator (fT3j) . However, the results shown in Figs. El - 
El illustrate the importance of a consistent treatment of pairing correlations by including an 
isospin invariant pairing interaction in the description of the ground state and the dynamics 
of lAR inN ^ Z nuclei. 



B. The Gamow- Teller resonance 



The Gamow- Teller resonance represents a coherent superposition of high-lying J'^ = 
1^ proton-particle - neutron- hole configurations of maximum collectivity associated with 
charge-exchange excitations of neutrons from orbitals with j = I + \ into proton orbitals 
with j = / — |. The GT operator reads 

T,^J = ESr±. (17) 

i=l 

The calculated GT strength distributions for ^^Ca, ^°Zr and ^osp]^ ^^^^ shown in Fig. (HI In 
addition to the high-energy GT resonance - a collective superposition of direct spin-flip 
{j = l + \^j = l — \) transitions, the response functions display a concentration of 
strength in the low-energy tail. The transitions in the low-energy region correspond to core- 
polarization (j = / ± I i = / ± |), and back spin-flip {i = l — \-^i = l + \) neutron- hole 
- proton-particle excitations. 

As in the calculation of the lAR, the effective RMF interaction is DD-MEl. In addi- 
tion, the standard values have been used for the parameters of the pion-nucleon interac- 
tion Lagrangian: = 138 MeV and /^/47r = 0.08, and the parameter of the zero-range 
Landau-Migdal force g' = 0.55 has been adjusted to reproduce the excitation energy of the 
GT resonance in ^''^Pb. The calculated resonances are shown in comparison with experi- 



mental data 

208 



391, 



41 



thick arrows) for the GTR excitation energies in 4«Ca |37|, 90Zr HQ, 

and 



42j . Although the residual interaction has been adjusted to reproduce the 
GTR excitation energy in ^°^Pb, we notice a very good agreement with experiment also 
for ^^Ca and ^^Zt. The adjusted value of g' in general depends on the choice of the RMF 
effective interaction. In the RRPA calculation with the NLl 
Ref. the GTR excitation energies of ^^Ca, ^°Zr and ^osp]^ were best reproduced by 
g' = 0.7, whereas approximately the same quality of agreement with experimental data was 
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43| effective interaction of 



obtained in Ref. [30| by using the NL3 jSjJ effective interaction and g' = 0.6. By using a 
set of density-dependent RMF effective interactions from Ref. we fiave verified tfiat, in 
general, there is a correlation between the value of the nuclear asymmetry energy at sat- 
uration 04, and the value of g' adjusted to reproduce the GTR. In order to reproduce the 
experimental excitation energy of the GTR, effective interactions with higher values of 04 
generally require higher values of g'. 



Gamow-Tel 
attention 29 , 



er states calculated in relativistic models have recently attracted considerable 
145^ . In particular, it has been shown that in a relativistic RPA 



3, y, 



calculation the total GT strength in the nucleon sector is reduced by 



2% in nuclear 



matter, and by ~ 6% in finite nuclei when compared to the Ikeda sum rule 

{S^T-S^^)=3iN-Z), (18) 

where S^± denotes the total sum of Gamow- Teller strength for the transition. The 
reduction has been attributed to the effect of Dirac sea negative-energy states, i.e. the 
missing part of the sum rule is taken by particle-hole excitations formed from ground-state 
configurations of occupied states in the Fermi sea and empty negative-energy states in the 
Dirac sea. The effect is illustrated in Fig. [71 where we display the running sum of GTR 
strength for ^°^Pb. The horizontal dotted line corresponds to the value 3(A^ — Z) = 132 
of the Ikeda sum rule. The solid and dashed lines corresponds the values of the GTR 
sum calculated from —00 to the excitation energy denoted on the abscissa. The big jump 
in the calculated GTR sum occurs, of course, when the main GTR peak at ~ 19 MeV 
is included. The RRPA calculation represented by the dashed line includes only positive 
energy ph configurations. Even extending the sum up to 60 MeV, the total sum amounts 
only to ^ 122, that is 8% less than the Ikeda sum rule. The total sum rule 3(A^ — Z) 
is exhausted by the calculated GT strength only when the relativistic RPA/QRPA space 
contains ph excitations formed from ground-state configurations of the fully or partially 
occupied states of positive energy, and the empty negative-energy states from the Dirac sea 
(solid line in Fig. [7j). The corresponding discrete GT spectrum is shown in Fig. |H| Two 
regions of excitation energies are shown. The panel on the right-hand side contains the 
positive energy irp — uh strength, with a pronounced Gamow- Teller peak at ^ 19 MeV. The 
left panel displays the negative energy spectrum built from na — vh transitions (a denotes 
a negative energy state). Even though these transitions are much weaker than the GTR, 
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there are many of them and their overall sum represents the strength missing in the positive 
energy sector. Similar results for the GTR sum rule are also obtained for ^^Ca and ^°Zr. 

For the case of open-shell nuclei we will compare the PN-RQRPA results with experi- 
mental data on Gamow- Teller resonances obtained from Sn(^He,t)Sb charge-exchange re- 
actions [40j. The GT strength distribution in ^^^Sn is shown in Fig. El Direct spin-flip 
transitions {uj = / + | ^ vrj = / — |) dominate the high-energy region above 10 MeV. The 
low-energy tail of the strength distribution corresponds to core-polarization (z/j = / ± | — »• 
vrj = / ± |), and back spin-flip (z/j = I — \ ~* ttj = / -|- |) transitions. 

The solid curve in Fig.lHlhas been calculated without including the T = proton- neutron 
pairing in the RQRPA residual interaction. The resulting high-energy GT strength dis- 
plays pronounced fragmentation. This fragmentation, however, is not induced by coupling 
to 2p — 2h configurations (not included in the present version of the RQRPA), rather it is 
caused by the splitting between different ph configurations. GTR configuration splitting (an 
appearance of two or more collective bumps with comparable intensities in the GTR strength 
function) was investigated in Ref. in the framework of the shell optical model. For Sn 
nuclei, in particular, this effect was predicted to occur as the valence neutron start to occupy 
the level with the highest j in the shell - /;.ii/2- The configuration splitting was attributed 
to the fact that the unperturbed energies of the {lgji2){)-g^i'^~'^ and (l/ig^g) (1^11/2)""'^ 
figurations are almost degenerate. The residual interaction removes this degeneracy and, as 
a result, the main GT component separates in two distinct peaks. The ground-state pairing 
correlations have a strong influence on the occupation of the l/iii/2 level, and therefore the 
energy spacing between the two peaks will depend on T = 1 pairing. For ^^'^Sn the calculated 
splitting of the GTR was 2.6 MeV |46|]. Subsequently, the fragmentation and splitting of the 
GTR in Sn nuclei was experimentally investigated in Ref. j40|. The theoretically predicted 
configuration splitting of the main GT component, however, could not be observed, since 
the total widths of the resonances of 5 — 6 MeV exceed the predicted splitting. The splitting 
of the main GT component shown in Fig. IHl (solid line) of ~ 3 MeV is in agreement with 
the result of Ref. In addition, we find a third direct spin-flip component in the region 
above 10 MeV, predominately based on the configuration (2(^3^2) (2(i5^2)^^- 

The other curves shown in Fig. El (dotted, dashed and dot-dashed) have been calculated 
with the inclusion of the T = proton-neutron pairing (fT^ in the RQRPA residual inter- 
action. In Ref. the overall strength parameter Vq of the interaction was adjusted to 
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measured half-lives of neutron-rich nuclei in regions where r-process path comes closest to 
the valley of stability: Vo = 230 MeV near = 50 and Vq = 170 MeV near = 82. In 
the present analysis, however, /?-decay half-lives are not calculated and we could not adjust 
Vq in this way. The GT strength functions shown in Fig. El have been calculated for Vq = 
200, 250 and 300 MeV, respectively. As one would expect, the inclusion of T = pairing 
has a strong influence on the low-lying tail of the GT distribution, i.e. on the region that 
contributes to /5-decay. Here, however, we are concerned with the main GT component. An 
important effect that we notice is the disappearance of configuration splitting between the 
two high-energy peaks. This happens because the T = pairing interaction does not affect 
configurations based on the ( 1(79/2 ) orbital (fully occupied), whereas it lowers configurations 
based on (partially occupied). Our calculation therefore shows that the 

T = proton-neutron pairing strongly reduces the predicted configuration splitting of the 
main high-energy GT component. In addition to the main GTR at ~ 16 MeV, however, 
part of the strength associated with direct spin-flip transitions is concentrated at ~ 10 MeV. 
Another interesting result is that the centroid of the GT strength composed of direct spin- 
flip transitions (z/j = / + | — > ttj = / — |) practically does not depend on the strength of 
the T = proton-neutron pp residual interaction. This is shown in the insert of Fig. El 
where we plot the calculated centroid of the direct spin-flip GT strength, as function of the 
strength parameter Vq. Since we could not use experimental data on the GTR centroid to 
adjust the strength of the T = pairing channel, the following illustrative calculations have 
been performed with Vq = 250 MeV. 

In a recent PN-RQRPA calculation Q] we have shown that the isotopic dependence of 
the energy spacings between the GTR and lAR provides direct information on the evolution 
of neutron skin-thickness along the Sn isotopic chain. It has been suggested that the dif- 
ference between the radii of the neutron and proton density distributions along an isotopic 
chain, could be determined from the measurement of the excitation energies of the GTR 
relative to the lAR. The ealeu.at.o. of Ref. Q. howeve. was perfo.tned w.th the NL3 
relativistic mean-field effective interaction. Although NL3 has become a standard in RMF 
calculations, like many other non-linear meson-exchange effective interactions it has a high 
nuclear matter asymmetry energy, and consequently it predicts rather large values for neu- 
tron density distribution radii. Here we repeat the calculation using the DD-MEl effective 
interaction. Among other ground-state nuclear properties, the parameters of DD-MEl have 
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also been adjusted to experimental data on differences between radii of neutron and proton 
distributions. 

In Fig. Uniwe display the calculated differences between the centroids of the direct spin-flip 
GT strength and the respective isobaric analog resonances for the sequence of even-even Sn 
target nuclei. For A = 112 - 124 the results of RHB plus PN-RQRPA calculation (DD-MEl 
density-dependent effective interaction, Gogny T = 1 pairing, T = pairing interaction 
fll4|) with Vq = 250 MeV, the Landau- Migdal parameter g' = 0.55), are compared with 
experimental data [40]. The calculated energy spacings are in very good agreement with 
the experimental values. As it has been emphasized in Ref. [s^, the energy difference 
between the GTR and the IAS reflects the magnitude of the effective spin-orbit potential. 
Many relativistic mean-field calculations have shown that the magnitude of the spin-orbit 
potential is considerably reduced in neutron-rich nuclei 47|, and a corresponding increase 
of the neutron skin has been predicted. In Fig. ^2 the calculated and experimental energy 
spacings between the GTR and IAS are plotted as a function of the calculated differences 
between the rms radii of the neutron and proton density distributions of even-even Sn 
isotopes (upper panel). Notice the uniform dependence of the energy difference between 
the GTR and IAS on the size of the neutron-skin. This means that, in principle, the 
value of Tn — Tp can be directly determined from the theoretical curve for a given value of 
-Eqt ~ -2'iAs- In the lower panel the calculated differences between neutron and proton rms 
radii are compared with available experimental data [4^ . The agreement between theoretical 
and experimental values suggests that the neutron-skin thickness can be determined from 
the measurement of the excitation energies of the GTR relative to IAS. 



IV. CONCLUDING REMARKS AND OUTLOOK 



In this second part of the work on the relativistic quasiparticle random phase approxi- 
mation (RQRPA), we have formulated the proton-neutron RQRPA in the canonical single- 
nucleon basis of the relativistic Hartree-Bogoliubov (RHB) model, for an effective Lagrangian 
characterized by density-dependent meson-nucleon couplings. The basic advantage of work- 
ing in the canonical basis is that, since it diagonalizes the density matrix, the canonical basis 
is always localized and so it can be used both for bound states and for the positive energy 
single-particle continuum. This is especially important in the description of low-lying excited 
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states and giant resonances in weakly bound nuclei far from stability. Prom the practical 
point of view, it is also considerably simpler to solve the RQRPA equations in the canonical 
basis, in which the relativistic Hartree-Bogoliubov wave functions can be expressed in the 
form of BCS-like wave functions. 

The inclusion of relativistic effective interactions with explicit density dependence of the 
meson-nucleon couplings provides an improved description of asymmetric nuclear matter, 
nuclear ground-states and properties of excited states, especially in nuclei far from stability. 
The PN-RQRPA model includes both the T — 1 and T — pairing channels, and the 
model space contains pair configurations formed from the fully or partially occupied states 
of positive energy in the Fermi sea, and the empty negative- energy states from the Dirac sea. 
It has been shown that the inclusion of configurations built from occupied positive-energy 
states and empty negative-energy states is essential in order to satisfy model independent 
sum rules. 

The proton-neutron relativistic QRPA provides a natural framework for the description 
of spin and isospin excitations. In this work the model has been applied to the analysis 
of charge-exchange modes: isobaric analog resonances and Gamow- Teller resonances. A 
very good agreement with experimental data has been obtained, both for doubly-closed 
shell nuclei and for open shell nuclei. The importance of a consistent treatment of pairing 
correlations has been demonstrated in the example of lAR and GTR in Sn isotopes. 

In addition to the lA and GT 0?iw excitations, the (p,n) spectra show clear evidence for 
the collective spin-flip dipole {L — 1,S — 1) Ihcu excitation and spin-flip quadrupole (L = 
2, 5" = 1) 2?iuj excitation. The broad (> 10 MeV) spin-flip dipole and quadrupole resonances 
can be interpreted as superpositions of three collective modes with spin-parity = 0~ , 1~, 
2~, and J"" = 1"*", 2+, 3"*", respectively. It has not yet been possible to experimentally resolve 
these resonances into the different angular momentum components. Recently it has been 
suggested that experimental data on charge-exchange spin-dipole resonances can be used to 
extract information about neutron-skin thickness in neutron-rich nuclei. 

The PN-RQRPA will be applied to the study of spin-multipole charge-exchange modes. 
An interesting question is whether such a self-consistent approach, which starts from the 
description of the nuclear ground-state, can simultaneously reproduce the dynamics of these 
high-lying and, in many cases, only weakly collective states. Another important topic will 
be the calculation of /3-decay rates in semi-magic neutron-rich nuclei. These nuclei are 
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particularly important for r-process nucleosynthesys. Their lifetimes are rather long and 
they determine the abundances of stable nuclides. Very neutron-rich nuclei are also very 
difficult to reach experimentally, and at present the only information on /3-decay rates is 
provided by theoretical calculations. This also emphasizes the importance of the comparison 
of the present approach with results of non-relativistic QRPA calculations. 
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FIG. 1: PN-RQRPA J'" = 0+ strength distributions. The excitations of the isobaric analog 
resonances are compared with experimental data (thick arrows) for ^^Ca js^], ^°Zr 
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FIG. 2: RHB plus proton-neutron RQRPA results for the isobaric analog resonances of the se- 
quence of even-even Sn target nuclei. The experimental data are from Ref. 4C|. 



FIG. 4: Same as in Fig. EJ for even-even Sn target nuclei with A = 120 — 132. 

FIG. 5: PN-RQRPA J"" = 0+ discrete spectra for ^^^Sn, ^^^Sn and ^^OSn. The panels on the left 
correspond to calculations performed without including pairing correlations. T = 1 Gogny pairing 
is included in the RHB calculations of ground states in the middle panels. The panels on the 
right-hand side display the results obtained by fully-consistent PN-RQRPA calculations, including 
the T = 1 proton-neutron pairing channel. 



FIG. 6: Gamow-Teller strength distributions for "^^Ca, ^^Zr and ^ospb. PN-RQRPA results are 



shown in comparison with experimental 
48Ca y, 90Zr 



mm,, and ^'J^Pb 



3S 



data (thick arrows) for the GTR excitation energies in 



42|. 



FIG. 7: The running sum of the GTR strength for Pb. The dashed line corresponds to a PN- 
RRPA calculation with only positive-energy ph configurations. For the calculation denoted by the 
solid line the RRPA space contains configurations formed from occupied states in the Fermi sea 
and empty negative-energy states in the Dirac sea. The total sum of the GT strength is compared 
to the model independent Ikeda sum rule (dotted line). 

FIG. 8: The PN-RQRPA strength distribution of discrete J"" = 1+ GT" states. The panel on the 
right-hand side contains the positive energy vrp — vh strength. The left panel displays the negative 
energy spectrum built from ttq — uh transitions. 



FIG. 3: PN-RQRPA J"^ = 0"*" strength distributions for even-even Sn target nuclei with A = 
108 — 118. For description see text. 
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FIG. 9: The Gamow- Teller strength distribution in Sn, calculated for different values of the 
strength parameter of the T = pairing interaction (|14() . The insert displays the corresponding 
calculated centroids of the direct spin- flip GT strenght. 



FIG. 10: RHB plus proton-neutron RQRPA results for the energy spacings between the Gamow- 
Teller resonances and the respective isobaric analog resonances for the sequence of even-even 
ii2-i24gj^ target nuclei. The experimental data are from Ref. 



FIG. 11: The proton-neutron RQRPA and experimental |4(j| differences between the excitation 
energies of the GTR and lAR, as a function of the calculated differences between the rms radii of 
the neutron and proton density distributions of even-even Sn isotopes (upper panel) . In the lower 
panel the calculated differences r„ — are compared with experimental data |43] • 
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